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Abstract 

Expressions are not functions. Confusing the two 
concepts or failing to define the function that is com- 
puted by an expression weakens the rigour of interval 
arithmetic. We give such a definition and continue 
with the required re-statements and proofs of the fun- 
damental theorems of interval arithmetic and interval 
analysis. 

1 Introduction 

Make things as simple as possible, 
but not simpler. 

Albert Einstein. 

The raison d'etre of interval arithmetic is rigour. 
Yet it appears that the most fundamental fact, some- 
times referred to as the "Fundamental Theorem of In- 
terval Arithmetic", is not rigorously established. The 
fact in question can be described as follows. 

Let e be an expression with (xi, . . . ,x n ) as an or- 
dered set of variables (i.e. a finite sequence of distinct 
variables). Let / be the function in lZ n — > 1Z that is 
computed by e. Let the result of evaluating e with 
intervals Ii, . . . ,I n substituted for X\, . . . , x n be an in- 
terval Y. Then 
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Although this fact is fundamental to everything that 
is done in interval arithmetic, we have failed to find 
in the literature a definition of what it means for an 
expression to compute a function. In Section 11.21 we 
review the literature that we consulted. 

In (Tl]) the interval Y is typically considerably wider 
than the range of function values. Interval analysis 
relies on the fact that, as I\,...,I n become narrower, 
the sides in (fTJ) become closer to each other. A theorem 
to this effect, such as 2.1.1 or 2.1.3 in [9] deserves to 
be called Fundamental Theorem of Interval Analysis 
rather than interval arithmetic. 

Both theorems should rest on the foundation pro- 
vided by a definition of the function computed by an 



expression. We give such a definition for sets; as in- 
tervals are sets, the definition applies to intervals as a 
special case. 

1.1 Expressions and functions 

An expression is an entity consisting of symbols; it 
is an element of a formal language in the sense of com- 
puter science. Some of these symbols denote opera- 
tions; others are constants or variables and denote reals 
or intervals, according to the chosen interpretation. 

Unlike an expression, a numerical function is an el- 
ement of the function space lZ n — > 1Z, for a suitable 
positive integer n. Variables only occur in expressions, 
where they can re-occur any number of times. Vari- 
ables do not occur in functions; in fact, the notion of 
"occurs in" is not applicable to functions in lZ n — > 7Z. 
Instead, a function in lZ n — > 1Z is a map from n-tuples 
of reals to reals; the elements of the n-tuples are prop- 
erly called arguments, rather than "variables" . 

An additional reminder of the need to distinguish 
between expressions and functions is that different ex- 
pressions can compute the same function. Yet another 
reminder is that there exist functions that are not com- 
putable, whereas all expressions are, like programs, in- 
structions for computations. Contrary to programs in 
general, expressions of the type of interest to interval 
arithmetic can be evaluated in finite time. Hence the 
functions computed by these expressions belong to the 
computable subset of functions. 

Of course, "the set of expressions" could be made 
precise by means of a formal grammar. For the purpose 
of this paper, it is sufficient to define an expression as 
follows. 

1. A variable is an expression. 

2. If E is an expression and if tp is a unary operation 
symbol, then ipE is an expression. 

3. If Ei and E% are expressions and if is a binary 
operation symbol, then E\ o Ei is an expression. 

To make the definition formal, we would have to spell 
out the appearance of the variables and of the operation 
symbols. 
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In whatever way expressions are denned, the result- 
ing set is disjoint from the set 7Z n — > TZ, whatever n 
is. What is needed to turn the Fundamental Fact ([TJ) 
into a theorem is to define "function computed by an 
expression" as a mapping from the set of expressions 
in n variables to 7Z n — > TZ. As observed above, this 
mapping is neither injective nor surjective. This map- 
ping can be called the semantics of the language of 
expressions. 

1.2 Remarks on the literature 

Moore [8] avoids the problem of defining the func- 
tion defined by an expression by not making the dis- 
tinction. As explained in the previous section, this is 
not correct. Jaulin et. al. [BJ, Theorem 2.2, assume 
that the problem is taken care of by composition of 
functions, but make unjustified simplifications. Com- 
position is indeed a promising approach, which we will 
pursue in Section [4] 

Neumaier [9] does distinguish between expressions 
and functions, but the expressions as he defined them 
fail to be computable. In fact, following the definition 
he gave in page 13, every real number is an element 
of the set of arithmetical expressions. The simplicity 
arises from the fact that all real numbers are defined as 
(sub)expressions. This introduces infinite expressions: 
whatever notation is chosen for the reals, most (in the 
sense of a subset of measure one) are infinite. In this 
way effective computability is lost. 

Moreover, Neumaier starts with an arithmetic ex- 
pression /, and then defines the interval evaluation 
of /, which he denotes by the same symbol /. To 
deal with partial functions, he introduced a NaN sym- 
bol, and the results of operations on this symbol. He 
then defined the restriction of / to its real domain 
D f = {x G Tl n | f(x) ^ NaN} to be the real eval- 
uation of /. We do not see the need for this indirect 
approach: partial functions are a perfectly natural and 
hardly novel generalization of functions that are total. 

Ratschek and Rokne also distinguish expressions 
from functions. In [T^ they refer to their earlier book 
[11) for a definition. This is a mistake, because on 
page 23, after a heuristic discussion of the connection 
between expression and functions, they refer to texts in 
logic and universal algebra for a definition. However, 
these assume that all functions are total. This is not 
always the case for the expressions of interest to inter- 
val arithmetic; consider for example \Jx. As only a few 
exotic varieties of logic allow function symbols to be 
interpreted by partial functions, it is better for interval 
arithmetic to use set theory as basis for its fundamen- 
tal theorems. In fact, these exotic varieties are subject 
to considerable controversy [2 [10] , so not suitable as a 



fundament for interval analysis. 

2 Set theory preliminaries 

This section establishes the concepts, terminology 
and notation for this paper. It is necessary because the 
present investigation is unusual in that all functions are 
what are usually called "partial functions" . To avoid 
having to qualify with "partial" every time a function 
is mentioned, we define "function" to mean what is 
usually referred to as "partial function". In other re- 
spects, we adhere closely to standard expositions of set 
theory, such as [31 Q] and standard introductions such 
as found in authoritative texts such as [7]. 

Definition 1 A function / consists of a source, a tar- 
get, and a map. The source and target are sets. The 
map associates each element of a subset of the source 
with a unique element of the target. 

The set of functions with source S and target T is 
denoted by the term S — > T. If a function f £ S —> T 
associates x G S with y G T, then one may write y as 
/(x). When only the association under / between x 
and y is relevant, we write x <— > y. 

Example 1 The square root is a function in TZ — > TZ 
that does not associate any real with any negative real 
and associates with x G 1Z the unique non-negative y G 
TZ such that y 2 — x if x > 0. 

The term /(x) is undefined if there is no y G T 
associated with x€S by feS^T. We take {f(x) | 
x G S} to mean 

{y G T | 3x G S and / associates y with x}. 

That is, {f(x) \ x G S} is defined even though f(x) 
may not be defined for every x G S. 

Example 2 {s/x \ x G TZ} is defined and is the set of 
non-negative reals. 

{x/y | x G {1} and y G TZ} is defined and is7Z\{0}. 

The subset of S consisting of x with which / G S — ► 
T associates a y G T is called the domain of /, denoted 
dom(/). If dom(/) = S, then / is said to be a total 
function, {/(x) | x G S} is called the range of /. 
We introduced the unusual terms "source" for S and 
"target" for T because of the need to distinguish them 
from "domain" and "range" . 

Definition 2 The set of functions with source S and 
target T is denoted S — > T and is called a "type" or 
"function space". 

Again, this differs from the usual meaning of S — > T, 
where it only contains total functions. To say that / 
"is of type" S -> T means that / G S -> T. 



Definition 3 Let f G S -> T and g G T -> [/ . T/ie 
composition go/ o/ / and g is the function in S — > Z7 
smc/i i/iaf g°/ associates x G 5 im'f/i z £ U iff there 
exists a y G T suc/i t/ia£ / maps x to y and g maps y 
to z. 

This is the conventional definition of composition. 
It requires the target of / to be the same set as the 
source of g. Because of this requirement it is not clear 
what composition Jaulin et. al. have in mind in [5J, 
Theorem 2.2. 

It follows from Definition [3] that the domain of def- 
inition of / o g is a subset of that of /. 

Example 3 / ogoh has {0} as domain if f G 72 — * 72 
is such that it maps x i— ► y/x, g G 72 — ► 72 is such that 
it maps x i— ► — x, and /i € 72. — ► 72. is such that it maps 
x i ► |x|. In other words, y/(— |x|) is undefined for all 

x G 72. except when x = 0. 

Let f € S —> T. The elements of 5 1 are called "argu- 
ments" of /. Note that if a function associates an x in 
S 1 with a y in T, it only so associates a single element of 
5. In that respect, all functions are "single-argument" 
functions. But S and T may be any sets whatsoever. 
Suppose / G 72" — * 72. Now the single elements in 
the source of /, the arguments of /, are n-tuples of 
reals. Thus we interpret the usual /(xi,...,x„) as 
/((xi, . . . ,x„)). 

Definition 4 Let f x G Si — > T x and f 2 G S 2 — > 

T2. TTie Cartesian product 0/ /1 and fi, denoted 
fi x /2, is a function in Si x S2 — > Ti x T2 
having domain dom(fi) x dom(/2) = {(xi,X2) 

xi € dom(fi) and x 2 G dom(/2)}, and mapping every 
(xi,x 2 ) indom(fi) x dom(/ 2 ) to </i(xi), / 2 (x 2 )) . 

Definition 5 Lei f be a function in S — > T . Lei 77 &e 
a total function in V(S) — > V(T). F is a set extension 
0/ / #f {/(x) 1 £ 1} C /or aH s«6sets X 

0/ S 1 . TTie tofaZ function in V(S) — ► V(T) with map 
X 1 ► {/(x) I x G X} is a se£ extension and is called 
the canonical set extension 0//. M^e mZZ use /(7?) to 
denote {/(x) | x G Z?}. 

3. Intervals are sets — interval exten- 
sions are set extensions 

As we saw, partial functions have set extensions that 
are total. This is of particular interest in numerical 
computation, where some important functions, such as 
division and square root, are not everywhere defined. 

In some treatments of interval arithmetic this leads 
to the situation in which division by an interval con- 
taining zero is not defined. This is not necessary: if 



one regards an interval as a set and an interval exten- 
sion as a set extension, then the interval extension is 
everywhere defined. This is the approach taken in [5], 
which will be summarized here. 

A well-known fact is that the closed, connected sets 
of reals have one of the following forms: {x G 72 | x < 
6}, {x G 72 I a < x}, {x G 72 | a < x < b}, as well 
as 72 itself. Here a and b are reals. Note that the 
empty subset of 72 is an interval also, as no ordering is 
assumed between a and b. 

The closed, connected sets of reals are defined to be 
the real intervals. They are denoted [—00,6], [a, 00], 
[a, 6], and [—00,00]. These notations are just a short- 
hand for the above set expressions. They are not meant 
to suggest that, for example, — 00 G [— 00, 6] = {x G 72 | 
x < b}. This is not the case because [—00, b] is a set of 
reals and —00 is not a real. 

The floating-point numbers are a set consisting of 
a finite set of reals as well as —00 and 00. The real 
floating-point numbers are ordered as among the reals. 
The least (greatest) element in the ordering is —00 (00). 
The floating-point intervals are the subset of the real 
intervals where a bound, if it exists, is a floating-point 
number. We assume that there are at least two finite 
floating-point numbers. As a result, the empty subset 
of 72 is also a floating-point interval. 

The floating-point intervals have the property that 
for every set of reals there is a unique least floating- 
point interval that contains it. This property can be 
expressed by means of the function □ so that OS is 
the smallest floating-point interval containing S C 72.. 
Given a real number x, we denote by x~ the greatest 
floating-point number not greater than x, and by x + 
the least floating-point number not less than x. 

By themselves, set extensions are not enough to ob- 
tain interval extensions. They need to be used in con- 
junction with the function □, as in the following defi- 
nition of interval addition: 

X + Y = U{z G 72 I 3x G X,y G Y.x + y = z} (2) 

for all floating-point intervals X and Y. Compared 
with a definition such as 

[a,b} + [c,d} = [{a + c)-,(b + d) + ], (3) 

(which is equivalent for bounded intervals) has the 
advantage of being applicable to unbounded intervals 
without having to define arithmetic operations between 
real numbers and entities that are not real numbers. 
Moreover, ([2|) includes the required outward rounding. 

Similarly to ([2]) we have 
Definition 6 

def 

X + Y = n{z G 72 I 3x G X,y G Y.x + y = z) 



X-Y = a{z £ TZ | 3a; £ X,y £ Y.z + y = x} 
def 

X*Y = U{z £TZ\3x £ X,y £ Y.x * y = z} 
def 

X/Y = D{z £ TZ I 3x £ X,y £ Y.z*y = x} 
def 

^X = a{y £ TZ I 3x £ X.y 2 = x} 

Theorem 1 The functions defined in Definition [6| 
map floating-point intervals to floating-point intervals, 
are defined for all argument floating-point intervals, 
and are set extensions of the corresponding functions 
from reals to reals. 

This is a summary of several results in [5] . 
Definition 7 Let I be the set of intervals. F £ I n — > I 
is an interval extension of f £ TZ n — > TZ iff F is the 
restriction to domain I n C TZ n — > TZ of a set extension 
of f . F is the canonical interval extension of f is de- 
fined to be F(B) = {/(x) | x £ B} whenever this is an 
interval. 

4. Semantics of expressions via set the- 
ory 

As all but a few exotic varieties of logic restrict func- 
tions to be total, we develop the semantics of expres- 
sions on the basis of set theory, even though most treat- 
ments of set theory also restrict functions to be total. 
However, as we have seen, functions in the usual set 
theory are easily generalized so that totality is not as- 
sumed. Modifying logic so that function symbols can 
be interpreted by partial functions has graver repercus- 
sions unm. 

Suppose that the expression e has the form e\ + 
e 2 and that e\ computes f\ : 7l m — > TZ and that e 2 
computes / 2 : 7Z n — ► TZ. In such a situation, Jaulin 
et al. [6] (Theorem 2.2), suggest that the function / 
computed by e is the composition of +, f\ and fa. 

But such a composition is not possible, as the types 
do not match, as required in Definition[3] We can make 
a composition if we form the Cartesian product of f\ 
and / 2 and if we make additional assumptions about 
e\ and e%- To prepare these assumptions we need the 
following definition. 

Definition 8 Let {v\, . . . , v n } be the set of variables in 
expression e. The variable sequence of e is (v\, . . . , v n ) 
if the first occurrences of the variables in e are ordered 
according to this sequence. 

Consider the special case where m = n and where 
e\ and e 2 have the same variable sequence. Let 8 £ 
TZ n -> TZ n x TZ n with mapping 

(xi,...,x n ) i-> ((xx,...,x n ),{xi,...,x n )) 



As will be shown in Lemma [TJ the function computed 
by ei + e-2 is +o (f 1 x fz)a8. The types of 8, f\ x fi, and 
+ do match: they are, respectively, TZ' 1 — > (TZ n x TZ n ), 
{TZ n x TZ n ) -> TZ 2 , and TZ 2 -► TZ. Thus it is clear the 
composition is defined and that its type is TZ n — > TZ. 

But it is of course a very special case if e\ and e 2 
have the same variables in the same order of first oc- 
currence. To further illustrate what is needed to define 
a composition of +, ex, and e 2 , consider another special 
case: e\ and e 2 have no variables in common, and their 
variable sequences are . . . , v m ) and (w\, . . . , w n ), 
respectively. As will be shown in Lemma [1] the func- 
tion computed by e\ + e 2 is again + o (f 1 x / 2 ) ° 8, 
except that 8 is in fZ m+n — > TZ m x TZ n and has as map 

\X\ , . . . , x m , yi , . . . , y n ) i ^ {{x%,. . .,x m ), (yi, . . . ,y n )) 

Now the types of 8, fi x / 2 , and +, are, respectively, 
n rn+n yjm x jjn^ (ftm x ^n) ^ n 2 ^ and ^2 _^ ^ 

Thus it is clear that the composition is defined and that 
its type is TZ' m+n -> ^. 

Finally, an example where the subexpressions share 
some, but not all variables. Consider the example 
where e\ is x * y, e 2 is j/ * z, e is ei + e 2 , and 
<5 e TZ 3 — > (T?. 2 x TZ 2 ) is such that 5 maps as follows: 
(xi,x 2 ,a;3) i-> ((x 1 ,x2), (x 2 ,x 3 )) for all xi,x 2 ,a;3 S 72.. 
Now the functions /i and / 2 computed by ei and e 2 
are the same function in TZ 2 TZ: it has as map 
(s,t) i— * s * £ for all reals s and £. Yet the function 
computed by e\ + e 2 does not have as map s * t + s * t: 
it is a different function, which is, however, described 
by the same formula + ° (/i X / 2 ) o <5. 

These three examples suggest how to define in gen- 
eral, for any pair (ei, e 2 ) of expressions and any domain 
D of interpretation, a "distribution function" that rep- 
resents the pattern of co-occurrences of variables in ex 
and e 2 . 

Definition 9 Given expressions e\ and e 2 with vari- 
able sequences (vi, . . . , v rn ) and (wi, . . . , w n ) , respec- 
tively. Let D be a set of values suitable for substitution 
for the variables. Let . . . , i p } and let {ji, . . . , j q } be 
a partition in {1, ... ,n} such that {w^, . . . iWi p } occur 
in ei and {vjj 1 Wj q } do not occur in e J^|. 

The distribution function i5 for the pair (ei, e 2 ) and 
D is the function in D m+q — > D m x D n that has as 
map 

(xi, . . . , x m , y n ,..., y jq ) i-* . . . , x m >, (yi, . . . , y„)) 
for all xi,... ,x m ,yi, . . . , y„ in D. 

1 Hence the variable sequence of any expression of the form 
ei (operation symbol)e2 is (vi, . . . , v m ,Wj 1 , . . . , Wj q ). 



Definition 10 An interpretation for an expression 
consists of a set D (the domain of the interpretation) 
and a map M that maps every n-ary operation symbol 
in the expression to a function in D n — > D. 

A set extension I' of I is said to be continuous if 
every symbol p is mapped to a continuous set extension 
of M(jp). I' is said to be canonical if every n-ary oper- 
ation symbol p is mapped to a canonical set extension 
ofM(p). 

The distribution function specifies enough of the 
way variables are shared between two expressions to 
support the central definition of this paper: 

Definition 11 Let e be an expression and let I be an 

interpretation that maps each n-ary operation symbol 
in e to a function in D n — > D, for n 6 {1,2}. We 
define by recursion on the structure of e, distinguishing 
three cases. 

Suppose e is a variable. The function computed by 
e under I is the identity function on D. 

Suppose e is ipei where if is a unary operation sym- 
bol. The function computed by e under / is fofi, where 
f is the function in D —> D that is the result of map- 
ping by I of ip and where f\ is the function computed 
by ei under I . 

Suppose e has the form e% o ei, where o is a binary 
operation symbol. Suppose 5 is the distribution function 
for (ei, ei) and D. The function computed by e under 

1 is o (/i x f%) o 5, where O is the result of mapping 
by I ofo. 

The definition assumes that no constants occur in 
expressions. We can simulate a constant by replacing it 
by a new variable and substituting the constant for that 
variable. In this way the definition does not suffer a 
loss of generality for expressions consisting of variables, 
constants, unary operators, and binary operators. At 
the expense of cumbersome notation (or sophisticated 
methods to avoid this), the function S can be extended 
to cover n-ary operation symbols with n > 2. 

The definition should conform to our intuition about 
expression evaluation. Suppose that D is the set of in- 
tegers, that the functions computed by e% and e 2 yield 

2 and 3, respectively. Then the definition should en- 
sure that the function computed by e% + e 2 yields 5 
when the interpretation maps + to addition over the 
integers. The following lemma confirms this intuition 
in general for arbitrary binary operation symbols. 

Lemma 1 Let e\ o ei be the expression in Defi- 
nition \11\ Suppose that {a\,...,a m ) G D m is 
substituted for (xi,...,x m ) and that (bi,...,b n ) € 
D n is substituted for (j/i, ...,y n ) . Let (ci,...,c q ) 
be such that (ai, . . . , a m , ci, . . . , c g ) is substituted for 
(xi , . . . , x m , j/i, ... , yq) . 



It is the case that 

/((ai, . . . ,a m ,ci, . . .,Cg)) 

= fl({ai, ■ ■ ■ ,a m ))Of2((bi, ■ ■ ■,&„)): 

where f is the function computed by e% o ei according 
to Definition \ll\ . 

Proof: 

/((ai,...,a m ,ci,...,c g )) = 

(Oo {fx x f 2 )o6)((a 1 ,...,a m ,c 1 ,...,c q )) = 

(O ° (/i x /2))<5((ai, . . . , Om, ci, . . . , c q ))) = 

(0° (fi x f2))((ai,...,a m ),(h,...,b n }) = 

0((fixf 2 )((( ai ,...,a m ),(b u ...,b n ))) = 
fi({a 1: . . . ,a m )) 0>/ 2 ((&i,. . .,&„)). 

Lemma 2 Let I be an interpretation for expression e 
and let I' be a set extension of I . Let f (f ) be the 
function computed by e under the interpretation I (V). 
Then f is a set extension of f. 

Though a minor lemma in set theory, the special 
case where the domains of / and /' are the reals and 
intervals respectively, it plays the role of the Funda- 
mental Theorem of Interval Arithmetic!!. 

Proof: We proceed by induction on the depth of the 
expression. Suppose the lemma holds for all expres- 
sions of depth at most n — 1. Let n be such that at 
least one of e± and ei is of depth n — 1 and the other is 
of depth at most n — 1. Suppose / has domain D and 
map M. Let e be e\ o ei and suppose that M maps o 
to O. Let 5 be the distribution function of ei and ei in 
that order. Let f\ and fi be the functions computed 
by ei and ei, respectively, under /. Let f[ and f% be 
the functions computed by e± and ei, respectively, un- 
der I . This gives as induction assumption that f[ and 
f 2 are set extensions of f\ and fi. 

Let / and /' be the functions computed from 
ei « 62 under interpretations / and respec- 
tively. Let Ax,...,A m ,Bi,,..,B n be subsets of D 
containing the elements a%, . . . , a m , 61, . . . , b n . Let 
c\, . . . , c q be such that S maps (ai, . . . , a m , c%, . . . , c q ) 
to ((ai,...,a m ), 6„)). 

Supposing that O' is a set extension of O, we have 

/((Ol,...,£t m ,Ci,...,C g » = 

/i(ai,...,a m )0/ 2 (&i,...,6„) G 
f[(A 1 ,...,A m )0'f! 2 (B 1 ,...,B n ) = 
f {Ax, . . . , A m , C\, . . . , C q ), 



2 Except that the statement in ]i inadvertently states instead 
the definition of interval extension. 



which is the function computed by e under Both 
equalities are justified by Lemma [T] 
Theorem 2 Let e be an expression with a variable se- 
quence (x\, . . . ,x n ). Let I be an interpretation for e, 
and I' a canonical set extension of L. Let f (f ) be the 
function computed by e under the interpretation I (V '). 
If each variable Xi occurs only once in e, then f is the 
canonical set extension of f . 

Proof: Following the same steps and notation as in 
the previous proof, we have 

f (A±, . . . , A m ,Bi, . . . , B n ) = 

(by Lemma [lj 

f{(A 1 ,...,A m )0'&(B 1 ,...,B n ) = 
(by the induction assumption) 

f 1 (A u ...,A m )0'f 2 (B u ...,B n ) = 

(using that f[ and f 2 are canonical set extensions and 
that ei and e 2 have no variables in common) 

{yeD\ 3 Vl e f 1 (A 1 ,...,A m ), 
3y 2 G h{B 1 ,...,B n ).y = y 1 Oy 2 } = 
{y G D | 3ai G Ay,..., 

3bi G B u ...,3b n G B n . 
V = h(ai, ■ ■ -,a m )Of 2 (b 1 , . . .,b n ) = 

f(ai,...,a m ,bi,...,b„)} = 
f(Ai, . . . , Am, Bi, . .., B n ). 

5. Continuous set extensions 

A fundamental fact in interval analysis can be stated 
intuitively as 

We can get arbitrarily close to the range of 
the point evaluation of an expression e by 
computing the interval evaluation of e with 
a sufficiently narrow interval. 

So far we were only concerned with interval arithmetic. 
This fact, being a continuity property, gets us into the 
realm of analysis. So it is here that interval analysis 
begins. 

As the validity of the statement and proof of such a 
property depends on a rigorous definition of the func- 
tion computed by an expression, it is wise to revisit the 
concepts and the theorems. 

Definition 12 Let T be a family of sets of D. A se- 
quence S = (S n } n £_\f of subsets of D converges with 
respect to T if it is nested, belongs to J- ' , and satisfies 
DneAT^n = { a }> w here a is an element of D. We say 
that the singleton set {a} is the limit of S . 



Definition 13 Let F G V(S) -> V(T), and let T x and 

be two families of sets of S and T , respectively. 
Let A = (A n ) n £N be any convergent sequence w.r.t. 
T\ with limit {a}. F is continuous w.r.t. T\ and F 2 
in {a} iff (F(A n )) n ^_\f is a convergent sequence w.r.t. 

Continuity is a very strong requirement. This raises 
the concern that no interesting examples might exist. 
The next lemma shows that this concern is unnecessary. 

Definition 14 Let f G K n -> K, and let \\.\\ be the 
Euclidean norm on lZ n . The function f is Cauchy- 
continuous at c £ lZ n iff for every e > there exists a 
5 > such that \\x — c\\ < S and x G dom(f) imply that 
|/0r) -/(c) | <e. 

A sequence (x^ieAA with Xi £ lZ n for all i G N is 
Cauchy-convergent to £ G TV 1 iff for every e > there 
exists an n such that ||£ — Xi\\ < e for all i > n. 

Lemma 3 Let f G TV 1 —tlZbe Cauchy- continuous 
at every x G dom(f) and suppose f has a canonical 
interval extension F. Then F is continuous w.r.t. the 
family of boxes oflZ n , and the family of intervals oflZ. 

Proof: Suppose that x is an element of lZ n , and that 
{Bnjrn^N' is a sequence of boxes in X™ that converges 
to x w.r.t. the family of boxes of TZ n . To prove that 
F is continuous w.r.t. the family of boxes of 7?™, and 
the family of intervals of 1Z, we have to show that the 
sequence (.F(-B„))„ eJ vconverges w.r.t. the family of in- 
tervals of 1Z. It is clear that this sequence is nested 
and belongs to the family of intervals of 1Z. So, we 
only need to show that C\ ne _\f F(B n ) is a singleton. In 
fact, 

□ F(B n ) = {f{x)}. 
The following inclusion is obvious: {f(x)} C 

r\neAf F ( Bn )- Lct V be ai1 element OI r\n£Af F ( B n)- 

This implies that for every n G Af, there exists x n 
in B n such that f(x n ) — y- Because (B n ) ne j^ is a 
nested sequence of boxes that intersect in {x}, the se- 
quence (x n )„ 6 7v' Cauchy-converges to x. Since / is 
Cauchy-continuous at x, we have f(x) = y. Therefore, 
(\ ne j^ F{B n ) C {/(x)}, which proves the lemma. 

Lemma 4 Let f G S — > T, and let T\ and Ti be two 
families of sets of S and T , respectively. Let F be a 
continuous set extension of f w.r.t. T\ and Ti and let 
A = (An)n,£j\f be a convergent sequence w.r.t. T\ with 
limit {a}. Then F({a}) = {/(a)}. 

Proof: As F is continuous w.r.t. T\ and Ti, 
(F(A n )) ne N is a convergent sequence w.r.t. J- 2 with 
limit, say, {b}. As F is a set extension of / we have 



that {/(x) | x £ Ai} C F(Ai), for all i e TV. As a e A: 
for all i € TV, we have that /(a) € {/(a:) | £ € Ai} for 
all i G TV". Hence /(a) e DieAT^^) = W- So we 
must have f(a) = 6. 

We are interested in interval extensions that are not 
canonical, yet are continuous. 

Starting from a family T of sets of a set D, we can 
construct a family of sets T n of D n , for any natural 
number n, by taking all the Cartesian products of any 
n sets in T . So, for any natural number n, and for 
any function F E V(D n ) — > V(D), we can study the 
continuity of F w.r.t. T n that was constructed from T . 
In this way, we treat the continuity of F by referring 
to T instead of T ri . 

In what follows, we suppose that the family of sets T 
of the domain D of an interpretation is given, and that 
the continuity of a set extension of an n-ary operation is 
based on this family. So, we will not use "w.r.t." from 
now on. In the case where D is TZ, T is the family of 
intervals in TZ. 

Definition 15 Let I be an interpretation with domain 
D and map M . A set extension I' of I is said to be 
continuous if every symbol p is mapped to a continuous 
set extension of M(p). V is said to be a canonical 
interval extension of I iff every symbol p is mapped to 
a canonical interval extension of M(p). 

Theorem 3 Let e be an expression. Let I be an inter- 
pretation for e, and let I' be a continuous set extension 
of I . Let f (F) be the function computed by e under 
the interpretation I (I 1 ). Then F is a continuous set 
extension of f . 

Proof: From Lemma [2j the function F is a set ex- 
tension of /. So we only need to prove that F is contin- 
uous. To do so, we proceed by induction on the depth 
of the expression e. The theorem holds when e has no 
subexpressions, that is, when e is a variable. In that 
case / and F are the identity functions, independently 
of I and V. The identity function in V(D) ->■ V(D) is 
continuous. 

This takes care of the base of the inductive proof. 
Let the induction assumption be that the theorem 
holds for all expressions of depth at most d — 1. Let 
e be the expression e\ o e 2 , where one of the subex- 
pressions has depth d — 1 and the other has depth at 
most d — 1. Suppose that the interpretation I has 
domain D and maps o to O. Let the interpretation 
I' have V(D) as domain and let it map o to O', a 
continuous set extension of O. Let 5 be the distribu- 
tion function with D for e\ and e 2 in that order. Let 
c\, . . . , c q be such that S maps (oi, . . . , a m , ci, . . . , c q ) 
to ((ai,...,a m ),(bi,...,b n )). 



Let F, F%, and F 2 be the functions computed 
under I' by e, ei, and e-x, respectively. Suppose 
that (A\) ieAr ,...,{Ai n ) ieM - and {B\) ie x, . . . , -B^); eA r 
are sequences of subsets of D that converge respec- 
tively to {ax},...,{a m } and {bi}, . . . ,{b n }. Ac- 
cording to the induction assumption F\ and F 2 
are continuous set extensions. This implies that 
({Fi(A\, . . . , A^), F 2 {B{, . . . , ))) iG /v converges to 
{(/i(ai, . . . , a m ), /s(i>i, • ■ ■ , b„))}, by Lemma E3 

Let (C^igA/ - be any such that 5(C l ) — 
{{A\,...,A i m ),{Bl...,Bi)) and such that {C l ) ieH 
converges to {{{a x , . . . , a m ), (h, . ■ ■ , b n )} 

We show that F is a continuous set exten- 
sion of / by showing that (F(C l ))i e j^r converges 
to {/((ai, . . . , a m , ci, . . . , Cq))}. To do so, we need 
to show that the sequence (F(C l ))i G Af is nested, 
and that C\i £ j^F(Ci) is the right value, namely 
{f({ax,...,a m ,ci,...,c q ))}. 

F{C l+1 ) = 

(by Definition HU 

(0'o(Fi x F 2 )oS)(C l+1 ) - 
(by application of 6) 

(o'o(f 1 xf 2 ))(«ai+ i , . . . , 4«>, (b?\ . . . , Birm = 

(by Definition H} 

0'((F 1 «4+ 1 , . . . , F 2 ({B?\ . . . , BI+ 1 )))) c 

(by the induction assumption and continuity of O') 

0'm((A[,. . .,Al n )),F 2 ((B{, . . . , B* m )))) = 

which proves that (F(C l ))i e j^f is nested. As for the 
convergence to the right value, we observe the follow- 
ing: 

n = 

(by Definition [TTj) 

f](0'o(F 1 xF 2 )o6)(C i ) = 

(by application of 5) 
D (<>' o (F X x F 2 ))(((A\,. ..,Ai),(Bi,..., Bi))) = 

(by Definition 

P| O'mdAl, O), F 2 ((B[,. . . , Bi)))) = 



(by continuity of O') 

0'« f) FrdAl AD), f) F 2 ({Bl = 

(by the induction assumption) 

0'«{/ 1 «o 1) . . . , a m ))}, {f 2 ((h, b n ))})) = 

(by Lemma SI 

{/ 1 ((a 1 ,...,a m ))}0{/ 2 ((6 1 ,..., &„))} = 

(because / is the function computed by e\ o e^) 

{/«ai, . . . ,a m ,ci, . . . ,c q ))}, 

which shows that F = O' o (i*\ x F2) o 5 is a continuous 
set extension of /, the function computed by e. 

Corollary 1 Let f 6 lZ n 1Z be the function com- 
puted by an expression e under an interpretation I that 
assigns Cauchy- continuous functions to the operation 
symbols in e. Let F be the function computed by e un- 
der the canonical interval extension of L. Let {Af)i e ^f 
be nested boxes converging to {a}. Then {F{A i )) il zj^ is 
a sequence of nested intervals converging to {/(a)}. 

In interval analysis, this corollary plays the role of Fun- 
damental Theorem. 

Proof: Since the image of any box by a Cauchy- 
continuous function is an interval, the interval exten- 
sion associated with each operation symbol is canonical 
(every Cauchy-continuous function has a canonical in- 
terval extension) . Using Lemma[3J these interval exten- 
sions are continuous. By Theorem [3J F is continuous. 
By Definition [13l (F(Ai)) ie _\f converges to {/(a)}. 

6 Conclusions 

The fact that the result of an expression evaluation 
in intervals gives a result that contains the range of val- 
ues of the function computed by the expression cannot 
be a mathematical theorem without a mathematical 
definition of what it means for a function to be com- 
puted by an expression. In this paper we give such a 
definition and prove the theorem on the basis of it. 

Another fundamental assumption in the use of in- 
tervals is that, as we make the intervals in an interval 
evaluation of an expression narrower, the interval re- 
sult gets closer to the range of values of the function 
computed by the expression. We use our definition to 
prove a theorem to this effect. 

Our starting point in all this is that intervals are sets 
and that, therefore, interval extensions of functions are 
set extensions of functions. The latter concept is an old 



one in set theory and is more widely applicable. Our 
definition and two main theorems are stated in terms 
of sets, so apply to intervals as special cases. 

This is of course only of interest to those who be- 
lieve in sets as foundation of mathematics. A radically 
different approach to the fundamental theorems of in- 
terval analysis is found in Paul Taylor's work (see for 
example |13]h Here the starting point is topology, ax- 
iomatically founded rather than set-theoretically. 

If it seems that our proposed foundations for inter- 
val methods are overly complex in comparison with the 
way they are given in the literature, we are comforted 
by Einstein's dictum: Make things as simple as possi- 
ble, but not simpler. 
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